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We generalize the Blonder-Tinkham-Klapwijk theory considering non-diagonal boundary condi¬ 
tions in the Bogoliubov-de Gennes scattering problem, to describe anomalous conductance features 
often reported for normal-metal/superconductor contacts. We calculate the differential conductance 
spectra showing that conductance dips, not expected in the standard formulation, are explained in 
terms of phase 7r-shift, between the bulk and the interface order parameter, possibly induced by 
a localized magnetic moment. A discretized model is used to give quantitative evaluation of the 
physical conditions, namely the polarization and transparency of the interface, needed to realize the 
phase gradient. 

PACS numbers: 74.45.-l-c; 74.25.-q; 74.25.Fy; 72.10.-d 


I. INTRODUCTION 

Since the introduction of the point contact spec¬ 
troscopy technique^ in 70’s, in which a micro-constriction 
is created pressing a metallic tip onto a superconduct¬ 
ing sample, the study of normal-metal/superconductor 
(N/S) junctions has represented an important means for 
the comprehension of several physical phenomena at the 
interface. The BTK theory^, formulated by Blonder- 
Tinkham-Klapwijk few years after, furnished a powerful 
tool to describe N/S contacts with transparency ranging 
from metallic to tunneling regime, the interface barrier 
strength being modeled using a Dirac delta potential of 
arbitrary amplitude. The theory, formulated in terms of 
Bogoliubov-de Gennes (BdG) equations^, provides the 
transmission and reflection coefficients and it succeeds 
in explaining the conversion of a quasi-particle current 
into a supercurrent, due to the Andreev reflection^, al¬ 
lowing accurate prediction of the experimental results 
about differential conductance spectra, energy gap and 
excess current. More recently, some modified BTK mod¬ 
els have been proposed in order to take into account 
spin polarization^ - —, diffusive contacts^, superconduct¬ 
ing proximity effect^ and thin ferromagnetic layers at 
the interfac o ^ ' - -- . These formulations have been moti¬ 
vated by several experimental evidence s 5 1 7 1 1 ° - — report¬ 
ing unusual conductance features, namely conductance 
dips and anomalous values of the zero-bias conductance 
(ZBC), not expected in the standard BTK model. 

From a mathematical point of view, including a Dirac 
delta potential within the BdG formalism, as done in 
the BTK approach, is equivalent to impose matching 
conditions for the scattering wavefunctions diagonal in 
the particle-hole representation^. However, off-diagonal 
boundary conditions are also mathematically allowed by 
the BdG formalism and they could account for the ap¬ 
pearance of anomalous features in the conductance spec¬ 
tra for N/S junctions. 

In this paper we extend the BTK approach to include 
non-diagonal boundary conditions in the BdG scattering 


problem, introducing an interface potential that mixes 
electron and hole components of the BdG state. We show 
that this potential describes the proximity effect at the 
interface and it is responsible for the formation of con¬ 
ductance dips in the differential conductance spectra of 
N/S contacts. The comparison of the differential conduc¬ 
tance curves with the experimental data suggests that a 
phase 7r-shift between the bulk and the interface order 
parameter, probably due to the presence of a localized 
magnetic moment at the interface, is at the origin of the 
conductance dips. The latter hypothesis is carefully an¬ 
alyzed by using a discretized model of the N/S junction. 

The paper is organized as follows: in Sec. II we formu¬ 
late the continuous model of the N/S interface introduc¬ 
ing a particle-hole mixing term in the scattering poten¬ 
tial. Off-diagonal boundary conditions are derived and 
the scattering coefficients are analytically determined. In 
Sec. Ill we show the differential conductance curves for 
N/S contacts obtained using generalized boundary con¬ 
ditions. We compare temperature evolution of conduc¬ 
tance spectra with the existing theoretical models and 
available experimental data. Possible phase shift effects 
at the interface are discussed. In Sec. IV we use a dis¬ 
cretized model to analyze the phase shift formation and 
evaluate the necessary physical conditions to observe it. 
Conclusions are given in Sec. V. 


II. MODEL 

We consider a one-dimensional N/S junction described 
by the BdG equations 

[H + V (x)] ip (x) = Eip (x ), (1) 

which completely define the quasi-particle state %)> ( x ) = 
(iq- ( x ), u\. ( x ), Vf ( x ), v\, (x)Y having excitation energy 
E above the Fermi energy Ep. The Hamiltonian H, 



2 


which describes the bulk properties of the junction is 


H = 


H 0 

—A*(x)ia y 


A {x)ia y \ 

-H* 0 )’ 


( 2 ) 


with 


Ho = 



E f I, 


(3) 


where I represents the identity operator in the spin space 
and a y is the Pauli matrix. We assume that the Fermi 
energy Ep and the effective mass m in the normal side 
of the junction (x < 0) are equal to those in the su¬ 
perconductor (x > 0), while the superconducting order 
parameter is taken of the form A(x) = A0(x), where 9(x) 
is the Heaviside step function. Differently from the stan¬ 
dard BTK treatment, we model the potential barrier at 
the interface (x = 0) by a particle-hole mixing operator 


V (x) 


Dot iUia y e lip 
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where Uo indicates the usual BTK barrier strength, while 
the term U\ describes the interfacial electron-hole cou¬ 
pling strength. The off-diagonal components of V (x) de¬ 
scribe the presence of a weak superconducting interface^ 
of negligible thickness compared to the coherence length 
of the superconductor. The variable tp represents the 
phase difference between the interface and the bulk su¬ 
perconducting order parameter. Maintaining arbitrary 
values of <p, a Josephson current Ij(<p) oc sin(y>)is is ex¬ 
pected to flow through the interface. The free energy of 
the system is expected to be minimized when Josephson 
current vanishes, i.e. for ip = 0 or w, the value ip = 0 
being a free energy minimum of the N/S junction. On 
the other hand, the value tp = it can become an energy 
minimum if a magnetic moment is formed at the interface 
(e.g., transition metals easily oxidize producing localized 
magnetic states). Indeed, in the presence of a sufficiently 
strong magnetic moment, the interfacial phase can be 
modified from 0 to 7r, the sign change of the interfacial 
order parameter following a mechanism similar to the one 
described in Ref. 0- 

In the following, we calculate the differential conduc¬ 
tance of the N/S junction by considering the generalized 
boundary conditions of the scattering problem. The wave 
function of an electron with spin a = {t,4-} coming from 


the N-side of the junction is given by: 
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Here the coefficients rj’4 and rJ/4 correspond, respec¬ 
tively, to normal reflection and Andreev reflection, while 
hk = \j2m (Ef + E) and hq = \j2m (Ep — E) indicate 
the electron and hole wave vectors. 

In the superconducting region, we have 


^s{x) = t\ 
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where the coefficients fj), correspond to the 

transmission as electron-like and hole-like quasiparticle 

with wave vectors Hk± = yj 2 m ( Ef ± \[~E 2 — A 2 ), the 
BCS2£ coherence factors being 

■ (7) 

The coefficients in Eqs. © and © can be determined by 
using the generalized boundary conditions for the wave 
functions at the interface: 

r N ( 0) = Vts(0) (8) 

|x=0 ^x^n\x=Q 

The matching matrix 

/ 0 in piw \ 

A = k F Zo\ M + k F Z 1 ( ^ J (9) 

contains a diagonal term in the particle-hole representa¬ 
tion with the usual BTK parameter Zq = 2mUo/(h 2 kF ), 
and off-diagonal terms of strength Z\ = 2mU\/(h 2 kp)- 
The Eqs. ©-© provide the simplest particle-hole mix¬ 
ing boundary conditions mathematically allowed by the 
BdG formulation. 

Using the above boundary conditions on the wave func¬ 
tions and the Andreev approximation (k+ = fc_ = k = 
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q = fcp) ; we find the following expression for the scatter¬ 
ing coefficients assuming the injection of a spin-up elec¬ 
tron from the normal side (the result doesn’t depend on 
the spin of the incoming process) 

4 . _ _ 4uv-2ie- xv Z 1 {u 2 -v 2 ) _ 

rfl 4 u 2 + AiuvZi cos ip + {u 2 — v 2 )(Zq + Z 2 ) 

AiuvZi cos ip + ( u 2 — v 2 )[Zo(2i + Zq) + Z 2 )] 
Au 2 + AiuvZi cosip + (u 2 — v 2 )(Z$ + Z 2 ) 

1 - Au — 2i(uZo — vZie~ lv ) 

e Au 2 + AiuvZi cos ip + (u 2 — v 2 )(Zq + Z 2 ) 

2i(vZo — uZ\e~ vp ) 

^ h Au 2 + AiuvZi cos ip + (u 2 — v 2 )(Zq + Z 2 ) ’ 

while the absence of spin-flip scattering implies r£ = rj, = 
t\ = f/ = 0. Once the scattering coefficients are ob¬ 
tained, we can calculate the differential conductance by 
the formula^ 
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where A a = \r^\ 2 and B a = \r%\ 2 are the Andreev re¬ 
flection and normal reflection probabilities, respectively, 
f(E) is the Fermi-Dirac distribution, while the notation 
a indicates the spin polarization opposite to a. 


III. RESULTS 

We first study the finite temperature conductance 
spectra of the N/S junction emphasizing the effects of 
the barrier strengths Zq, Z i and of the phase ip. In Fig- 
ure[T]we show the normalized conductance G/G n vs e/A, 
with G n = G(eV A), for different values of Zq and 
Z\ at a fixed temperature T = 1.6 K, computed by using 
Eq. (fill) . Two cases are considered: (i) <p = 0, shown in 
the left panels; (ii) ip = tt, shown in the right panels. In 
each plot, different curves correspond to different Zq val¬ 
ues (ranging from 0 to 2), while Z\ is fixed as labelled. 
For p = 0 and Z\ = 0 (Figure 1(a)) the usual BTK 
behavior is recovered. In this case, the zero-bias conduc¬ 
tance is suppressed as Zq is increased, while two peaks 
at e/A fts ±1 appear. Fixing Z\ = 0.5 (Fig. 1(c)), we 
observe a reduction of the amplitude of the zero bias 
conductance feature compared to the Z i = 0 case. The 
difference between the conductance lowering induced by 
Zq and the peculiar amplitude reduction induced by Z\ 
is evident: while the increasing of Zq induces a zero-bias 
conductance minimum, a tendency to increase the zero- 
bias conductance is observed by rising Z\ to 1.0 (Fig. 
1(e)), Zi = 1.5 (Fig. 1(g)) and Z x = 2 (Fig. l(i)). 

A different scenario is observed for ip = tt: the effect of 
moderate values of Z\, namely Z\ = 0.5 (Fig. 1(d)) and 
Z\ = 1 (Fig. 1(f)), combine with Zq to give a relevant 
effective barrier strength leading to a strong suppression 
of the sub-gap conductance up to fully gapped spectra. 



FIG. 1: (Color online) Normalized differential conductance 
curves, G/G n vs e/A, calculated at T = 1.6K from Equation 
(11). The different curves are obtained for distinct values of 
Zq, Z\ and ip (values in the panels). 


For Zi = 1.5 (Fig. 1(h)) and Z\ = 2 (Fig. 1(1)) an evi¬ 
dent zero-bias peak surrounded by two dips at e/A « ±1 
appears. Such ZBC peak exists for all Zq values in the 
range [0,2], the junction transparency reduction having 
effect only on the peak amplitude. 

All the conductance structures presented above (com¬ 
ing from the generalized boundary conditions) cannot 
be recovered within the standard BTK approach (except 
for the case Z i=0). Moreover, the interface potential 
given in Equation (U) can be further generalized to in¬ 
clude, spin-orbit interaction in the plane perpendicular 
to the transport direction, local magnetic moments and 
triplet or non-centrosymmetric superconducting correla¬ 
tions. The above complications make the interface po¬ 
tential an off-diagonal differential operator of the form 
B(x, d X ’ V ’ Z )6(x) acting on the Nambu space which induces 
an extended class of particle-hole mixing boundary con¬ 
ditions. Extending the BTK theory along this direction 
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produces analytic results for the scattering coefficients 
which can be directly employed to explain anomalous 
conductance spectra. 

In Figure 2 we compare the temperature evolution of the 
conductance spectra as obtained for ip = ir in the gen¬ 
eralized BTK model introduced above, with the two-gap 
models. In particular, in Fig. 2(a) we show theoretical 
curves calculated in the temperature range between 0.1 K 
and 5.1 K by assuming Z o =0.35, Z i=2, Ajv&= 1.5 meV, 
ip = 7r, while in Fig. 2(b) the conductance curves are ob¬ 
tained by considering the parameters Z= 0.14, A^O.99 
meV, A2=0.47 meV, in the two-gap model. Both models 
can be used to reproduce (black solid lines in Fig. 2(a) 
and 2(b)) the experimental data reported for Cu/Nb con¬ 
tacts in Ref. F° r both the theoretical curves is nec¬ 
essary to fix a temperature value of 0.9A', that is lower 
than the bath temperature and it has been motivated as 
the result of non-equilibrium effects^ as well as in terms 
of other physical effects^. However, the temperature 
evolution of the conductance spectra shows appearance 
of different low temperature features at zero bias as well 
as at the gap edge. Moreover, the temperature depen¬ 
dence of the structures, namely ZBC and Edip , evolve 
differently. By rising the temperature, a non-monotonic 
evolution of Ed ip (Fig. 2(c)) and a faster reduction of 
ZBC (Fig. 2(d)) is observed for the case of Fig. 2(b). 
As a consequence, very low temperature experiments are 
necessary to distinguish the two models, and thus under¬ 
standing the physical origin of the anomalous conduc¬ 
tance features observed in the point contact experiments 
in N/S devices. It is worth to notice that both scenarios 
recalled by the two models are physically plausible. In¬ 
deed, Nb and Cu oxides are known to exhibit magnetic 
correlations that could realize effective local polarization 
enabling a phase shift at the interface (discussed below); 
on the other hand, formation of a (proximized) weak su¬ 
perconducting layer at the N/S interface is always possi¬ 
ble. 

Up to now we have assumed that the interface phase 
ip = 7r is the result of an emergent magnetic moment at 
the interface even though the Hamiltonian model (Eq. 
(3)) does not include explicitly magnetic correlations. In 
the following, we prove that also modifying the interface 
potential V (x) (Eq. fl3])) by adding a localized magnetic 
potential of the form U 2 a z ® a z S(x), with dimensionless 
strength Z 2 does not change the peculiar shape of the 
conductance given in Fig. [2] Assuming an incoming spin- 
up electron, the expressions of the Andreev and normal 
reflection become: 

I Auv — 2ie~ l<p Zi(u 2 — v 2 ) 

h Au 2 + AituvZi cos ip + (u 2 — v 2 )Z+ 

^ AiuvZi cos (p + ( u 2 — v 2 )[Z+ + 2i(Zo + Z 2 )} 

4 u 2 + AiuvZi cos p> + ( u 2 — v 2 )Z 2 

where Z 2 = Z 2 + (Zq + Z 2 ) 2 . Analogous expressions, 
characterized by the parameter Z 2 _ = Z 2 + ( Zq — Z 2 ) 2 , 
are found for r' h and r / in the case of the scattering 
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FIG. 2: Temperature evolution of the normalized conduc¬ 
tance curves obtained for (a) generalized BTK model with 
particle-hole mixing boundary conditions and for (b) two-gap 
mode£ for the parameters Z=0.14, Ai=0.99 meV (bulk gap), 
A2=0.47 meV (proximized gap). The thick black line in the 
two plots represent the best fit for Cu/Nb experimental data 
reported in Ref. Q], (c) Comparison of the temperature evo¬ 
lution of the energy position of the conductance dips, Edip vs 
T, obtained from (a) and (b). (d) Comparison of the temper¬ 
ature evolution of the ZBC, obtained from (a) and (b). 



e/A 


FIG. 3: Normalized differential conductance curves, G/G n 
vs e/A calculated by considering the inclusion of a localized 
magnetic contribution to the potential V ( x ). Parameters are 
T = 1.6 K for Zq = 0, Z\ = 2 and ip = n, with Z 2 6 [0, 0.5]. 


problem of a spin-down electron coming from the nor¬ 
mal side. It is also worth mentioning that the magnetic 
moment makes different the Andreev (and the normal) 
reflection coefficients of scattering processes originated 
by electron-like quasiparticles of opposite spin projection 
(i.e. r^j h ^ r \ih)i while spin-flipping reflection processes 
are not allowed by the Hamiltonian structure. In Fig. 
[3] the effect of Z 2 on the conductance curves is shown. 
Increasing Z 2 , the amplitude of the zero-bias peak is 
lowered, while the conductance features are almost the 
same for Z 2 £ [0, 0.5] except for a small reduction of the 
ZBC peak amplitude. This shows that adding a localized 
magnetic term does not change the conductance shape. 
The latter property derives from the fact that the mag¬ 
netic correlation only renormalizes the interface potential 
(Zq —> Zq ± Z 2 ) producing a spin-sensitive effective bar- 
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FIG. 4: Discretized model of the N/S junction consisting of 
N (odd) sites: (N+l)/2 normal sites (A = 0) and (N-l)/2 
superconducting sites (A ^ 0); magnetic (F) e non-magnetic 
((7) potentials are present at the interface site i = (N + l)/2. 
The hopping parameter t is homogeneous along the system. 


rier, whose effects are difficult to be distinguished from 
the ones expected in non-magnetic case. This explains 
why the presence of an emergent magnetic moment at 
the interface is difficult to be experimentally confirmed 
by means of point contact spectroscopy. 


IV. DISCRETIZED MODEL 


As we have seen in Figs. 1 and 2 the conductance dips 
appear in the generalized BTK approach for ip = tt and 
this phase value can be associated to a localized magnetic 
moment. In fact, the presence of a localized magnetic 
moment at the interface can make a phase gradient of 
7 r energetically favorable. In order to identify the phys¬ 
ical conditions (interface polarization and transparency) 
to realize the 7r-shift, we consider a discretized formula¬ 
tion that allows to describe spatial dependent potentials 
without increasing the computational complexity. We 
model a system with an odd number of sites N in which 
( N — l)/2 sites are used for both the normal and the su¬ 
perconducting side, while one normal site with magnetic 
(r) and non-magnetic (U) potentials is assigned to the 
interface (see Fig. 4). The nearest-neighbor hopping pa¬ 
rameter t = h 2 /(2ma 2 ), expressed in terms of the sites 
distance a, is assumed homogeneous and it used as energy 
unit, fixing t ~ 10-Ajvb in order to have £Nb — 10-a. Tem¬ 
perature is measured in dimensionless units r = ksT/t. 
Hereafter, we set t = 16.2 meV as the energy cut-off of 
the theory, the latter being of the same order of mag¬ 
nitude of the Debye energy hi 2/j. This choice guaran¬ 
tees that only states with phonon-mediated attraction 
(i.e. with e £ [0,Ml£>]) are retained, simultaneously en¬ 
suring the long wavelength limit of the considered wave 
functions. Under these assumptions, the relevant wave 
functions present a De Broglie wavelength greater than 
the lattice constant a, while the associated eigenvalues 
defines a near-parabolic energy dispersion. 

The discretized version of the BdG equations, in the 
presence of a Zeeman term r(x)<j z added to the single¬ 
particle Hamiltonian H$ to account for the effective po¬ 
larization at the interface, correspond to the following 
matrix equations (V i £ [1, N]): 

+ T (*$ + *&) = 0 , ( 13 ) 





i r 


FIG. 5: Numerical results obtained in the discretized model. 

(a) Temperature dependence of the bulk superconducting gap 
for different systems, namely N = 15, N = 51, N = 71, as¬ 
suming transparent barrier, U = 0. Numerical data, normal¬ 
ized to the low temperature value Ao, are compared to theo¬ 
retical behaviour expected in the BCS model for r c = 0.0485. 

(b) Spatial dependence of the superconducting gap for 71-sites 
and 51-sites systems, by assuming U — 0, a constant BCS cou¬ 
pling, A i = A, and dimensionless temperature r = 0.025. Dif¬ 
ferent curves correspond to different values of T. Lines refer 
to data obtained for 71-sites system; scattered symbols refer 
to 51-sites system, rescaled to compare the data sets. The 
arrow indicates the region where the superconducting gap is 
calculated self-consistently. (c) Spatial dependence of the po¬ 
larization calculated for N = 71 for different T values, with 
(7 = 0 and r = 0.025. (d) Effect of the barrier strength U : the 
superconducting gap vs T is evaluated at the site i = 25 (of 
51), at r = 0.025 for three different transparency conditions 
(U = 0, U = 0.5, (7=1), the magnetic site being located at 
i = 26. 


where (a = ±, ± =t / 4-) 
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( CT ) 


f Si-E + aTi a A, 

V cA* -Si-E + aTi 


, ( 14 ) 


while T = —t a z . Here Ej/i = 2 + U5i^+ 1)/2 is the 

energy of the *-th lattice site, while = (u a ,i, 
is the discretized BdG state in the absence of spin-flip 
scattering. Fj/t = F<5 Z ,(jv+i )/2 is the site dependent Zee- 
man energy that we take different from zero only at 
the interface site. Using Dirichlet boundary conditions 
= 0, we get electron-like eigenstates 
n t 

* (a ' n) = v S) ( 15 ) 

i= 1 

associated to positive energy eigenvalues (e„ > 0), with 
Aj = ..., Sjsr^y. The spatial dependence of the su¬ 

perconducting gap is computed as^i 
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the sum being calculated for e n G [0, fifi.o]. The at¬ 
tractive phonon-mediated local potential A^ is assumed 
constant^ (A* = A) also in the normal side of the junc¬ 
tion to take into account for the proximity effect and the 
interdiffusion of atoms belonging to the N- and S-side. 
We consider the bulk superconducting gap A bulk at the 
center of the S-region, in order to avoid finite size effects. 
A bulk is self-consistently computed using Eq. msD with 
accuracy better than 1% starting from A, = 0.15t in S, 
and by fixing A = 2.4 t. The polarizing effect of the mag¬ 
netic site at the interface can be quantified defining the 
site-dependent polarization Pi = (n^— rii^/frii^+ni^), 
where 


ni^ + rii t i — + Kui*| 2 (l ~ fn) +(17) 

n 

+ 51 [l«S| 2 /n + 1*45 l^ 1 _ /") 


and 


ni,t - n.,4 = 5Z [\ u u\ 2 fn - 1*451^ 1 - fn) + (18) 
n 

+ 54 [ - i u JJi 2 /" + 1*451 2 ^ - /«) 


with f n = f(e n ) a shortened notation standing for the 
Fermi-Dirac distribution. In order to capture the bulk¬ 
like behaviour using a finite size system, the system 
size N has been progressively increased from N = 15 
to A = 71, while monitoring the temperature depen¬ 
dence of A bulk- The results of this analysis are shown 
in Figure 5(a) where normalized values of A \, u ik are pre¬ 
sented as a function of the dimensionless temperature 
r. For a system size of N = 15 a size-induced suppres¬ 
sion of the superconducting gap is observed, this effect 
being more evident close to the transition temperature 
t c . Increasing the system size up to N = 51 produces 
a A bulk vs r curve very close to the one obtained for 
the N = 71 case, signaling that the bulk limit of the 
interface model has been reached. The temperature evo¬ 
lution of A bulk for N = 71 has been compared with the 
BCS behaviour giving a dimensionless critical temper¬ 
ature r c = 0.0485, corresponding to a niobium critical 
temperature Tff b ~ 9.1 A'. In Figure 5(b) we present the 
spatial dependence of the superconducting gap (for sys¬ 
tem size N = 51 and N = 71) fixing the Zeeman energy 
r of the magnetic potential in the range [—2.5,0], while 
taking U = 0 (transparent interface) and r = 0.025 (i.e. 
T ~ Tff b / 2). In order to compare spectra obtained for 
systems with different size, the data referring to N = 51 
have been rescaled. For —1.0 < T < 0 we observe or¬ 
dinary proximity effect where finite superconducting or¬ 
der parameter is induced in the N-side on a length of 
about 10a — f,Nb- For sufficiently strong magnetization 
r < —1.0 negative order parameter is induced on the 
same length scale. Reduction of the order parameter on 


the right border is due to the S/vacuum interface. In 
Figure 5(c) we show the spatial dependence of the po¬ 
larization, calculated for T G [—2.5,0] and U = 0. The 
polarizing effect of the localized magnetic moment asym¬ 
metrically extends on a distance of about 20a. In the 
superconducting side the induced polarization is inverted 
for large T values^ (r < —2.0). The general aspect of 
the polarization curves evidences Friedel density oscilla¬ 
tions. We have also verified the effect of barrier strength 
on the inversion of the superconducting order parame¬ 
ter. In Figure 5(d), we show the gap value calculated (at 
site i = 25) in proximity of the interface, for a system 
size N = 51, as a function of T, with enhanced reso¬ 
lution (step 0.1). For reduced transparency (U > 0) a 
larger magnetic moment is necessary to induce the in¬ 
version of the interface order parameter. The analysis of 
the pairing potential A* shows that, in presence of a local 
polarization at the interface, a phase gradient tp = it can 
be stabilized. For relatively transparent junctions (i.e. 
described by small values of U) the sign change of the 
interface order parameter can be obtained with moder¬ 
ate polarization strength, while strong polarization val¬ 
ues are needed for opaque interface with higher values of 
U. Thus the probability to observe an hidden magnetic 
moment at the interface accompanied by a phase gra¬ 
dient is enhanced in transparent systems. The physical 
origin of a local magnetic moment at the Cu/Nb inter¬ 
face probably resides in many-body effects which can be 
accounted for in the framework of the Anderson impurity 
model^i. 


V. CONCLUSIONS 


We have generalized the BTK theory to include 
particle-hole mixing boundary conditions in the scatter¬ 
ing problem, reporting analytic results for the scattering 
coefficients. We calculated the finite-temperature differ¬ 
ential conductance spectra for N/S junctions showing the 
formation of conductance dips in the case of a phase tt- 
shift at the interface. We demonstrated that the tem¬ 
perature evolution of the conductance spectra can dis¬ 
criminate the physical origin of the conductance dips, 
either the formation of a localized magnetic moment or 
the presence of a weak proximized superconducting layer 
at the interface. According to the analysis, a localized 
magnetic moment could make a sign change of the su¬ 
perconducting order parameter energetically favorable. 
Finally, we used a discretized model to determine the 
necessary physical conditions under which the 7r-shift is 
realized: transparent interfaces can easily sustain a phase 
gradient as the effect of a weak interface magnetization, 
while for reduced transparencies a relative strong local¬ 
ized magnetization would be necessary. 
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